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Abstract 

The  demonstrated  advantages  of  the  probabil¬ 
ity  density  function  (pdf)  method  over  conven¬ 
tional  turbulence  models  have  made  it  a  most 
promising  tool  for  computing  the  complex  flow 
fields,  involving  turbulence,  mixing,  and  chemical 
reaction,  encountered  in  engineering  applications, 
especially  practical  combustion  devices.  There  are 
three  aspects  to  the  present  study:  i)  swirling 
flows,  an  essential  ingredient  in  practical  combus¬ 
tion  devices,  have  not  been  studied  previously  with 
the  pdf  method.  In  the  present  study,  the  pdf 
method  is  extended  and  applied  to  a  coaxial 
swirling  jet  flow  and  the  results  are  compared  with 
detailed  experimental  data;  ii)  the  recently  devel¬ 
oped  stochastic  dissipation  model  (Pope  and  Chen, 
1990),  hitherto  applied  only  to  homogeneous  turbu¬ 
lence  and  self-similar  flows,  is  applied  to  and  vali¬ 
dated  for  a  general  developing  shear  flow  as  repre¬ 
sented  by  foe  coaxial  swirling  jet  flow  considered. 
The  results  are  also  compared  against  computa¬ 
tions  using  a  previous  model  based  on  the  mean 

dissipation  (Anand,  Pope,  and  Mongia  1989);  iii) 

detailed  velocity  statistics,  conditional  upon  foe  jet 
from  which  the  fluid  originated,  are  calculated  and 
compared  against  measurements,  demonstrating 
that,  in  addition  to  other  advantages,  foe  pdf 
method  can  provide  considerably  more  description 
of  foe  flow  than  conventional  models,  and  yet  is 
computationally  tractable.  The  results  are  in  very 
good  agreement  with  data. 


*  Member,  AIAA 
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Introduction 

Several  previous  studies,  most  of  them  recently 
reviewed  by  Pope,1  have  established  foe  probabil¬ 
ity  density  function  (pdf)  method  as  a  prime  can¬ 
didate  for  computing  complex  turbulent  flow  fields 
in  which  mixing  and  reaction  are  occurring.  Unlike 
foe  currently  used  turbulence  models,  the  most 
important  processes,  namely  convection  (by  both 
mean  and  fluctuating  velocities)  and  reaction,  ap¬ 
pear  in  closed  form  in  the  joint  velocity-composi¬ 
tion  pdf  method  and  need  not  be  modeled.2  The 
conventional  turbulence  models  such  as  k-e  and 
Reynolds-stress  models  (RSMs),  based  on 
Reynolds-averaged  equations,  use  gradient  diffu¬ 
sion  models  for  turbulent  transport  (convection  by 
fluctuating  velocities).  More  importantly,  foe 
models  are  incapable  of  treating  finite-rate  multi- 
step  reactions.  The  ability  to  treat  these  reactions  is 
essential  for  computing,  understanding,  and  con¬ 
trolling  pollutant  emissions  that  has  become  ever¬ 
more  important  for  current  and  future  land-based 
and  aerospace  applications.  The  present  study  is 
another  step  in  foe  development  of  foe  pdf  method 
as  a  design  and  analysis  tool  for  practical  systems. 

The  pdf  method  has  been  applied  to  a  number 
of  reacting  and  nonreacting  flows  including  two- 
dimensional  recirculating  flows^  and,  in  conjunc¬ 
tion  with  a  Reynolds-averaged  finite-volume 
method,  to  two-  and  three-dimensional  time-de- 
pendent  flows.5'6  Swirling  flows  (flows  with  ap¬ 
preciable  bulk  swirl)  are  an  essential  feature  in 
many  practical  combustors  such  as  gas  turbine 
combustors.  Swirling  flows  have  not  been  studied 
previously  with  foe  pdf  method,  although  there  are 
no  theoretical  restrictions  for  applying  the  pdf 
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method  for  such  flows.  In  the  present  study,  the 
necessary  extensions  to  the  solution  algorithm  are 
made  to  include  the  treatment  of  swirl.  For  the 
boundary-layer  type  flow  under  consideration,  a 
salient  feature  of  this  extension  is  the  inclusion  of 
the  mean  axial  pressure  gradient  in  the  calcula¬ 
tions. 

It  is  well  known  that  the  k-e  model,  which  is 
the  most  widely  used  turbulence  model,  is  unsatis¬ 
factory  for  swirling  flows7  due  to  effects  of  stream¬ 
line  curvature,  although  the  model  performs  rea¬ 
sonably  well  for  nonswirling  jet  flows.  Several 
modifications  to  the  k-e  model  to  account  for  cur¬ 
vature  effects,  through  corrections  based  on  the 
Richardson  number,  have  been  investigated  and 
reported  in  die  literature.  The  corrections  are  typi¬ 
cally  made  through  the  equation  for  mean  dissipa¬ 
tion  as  a  means  of  correcting  the  turbulent  trans¬ 
port.  However,  these  modifications  have  not 
proved  to  be  uniformly  successful  over  a  range  of 
swirling  flows.  In  this  context,  it  is  necessary  to 
verify  the  performance  of  the  pdf  method  for 
swirling  flows. 

A  serious  deficiency  in  the  joint  velocity-com¬ 
position  pdf  method  is  that  the  joint  pdf  contains 
no  information  on  the  time  or  length  scales  of  the 
turbulent  motions.  However,  given  the  time  scale 
information,  the  method  provides  a  complete  one- 
point  description  of  the  joint  velocity-composition 
statistics.  For  several  shear  flows,  previous  stud¬ 
ies®"11  have  assumed  foe  turbulent  time  scale  to  be 
uniform  across  foe  flow.  The  turbulent  time  scale 
is  defined  to  be  foe  ratio  of  the  turbulent  kinetic 
energy  to  foe  mean  dissipation  rate  (t=k/<e>).  The 
assumption  of  a  uniform  time  scale  is  not  valid  in 
foe  developing  regions  of  foe  shear  flow.12  The  al¬ 
ternatives  have  been  to  solve3"6  foe  standard 
model  equation  for  <e>,13'14  or  to  solve  a  model 
equation12'13  for  foe  turbulent  frequency  <co>  = 
1/t. 

The  model  developed  recently  by  Pope  and 
Chen16  eliminates  the  need  to  supply  the  time  scale 
information  to  foe  pdf  equation.  This  is  done  by 
considering  the  joint  pdf  of  velocity,  composition, 
and  foe  instantaneous  dissipation  e.  In  addition  to 
effecting  complete  closure  of  foe  joint  pdf  equation, 
foe  inclusion  of  e  in  the  joint  pdf  leads  to  more  real¬ 
istic  modeling.  For  example,  foe  model  allows  for 
multiple  time  scales  rather  than  a  single  mean  time 
scale,  accounts  for  internal  intermittency,  has  po¬ 
tential  to  account  for  large  scale  structures,  and  bet¬ 
ter  accounts  for  foe  influence  of  the  origin  and  his¬ 
tory  of  the  fluid  particles  on  foe  local  turbulent 
structure.  The  model  was  constructed  with  refer¬ 


ence  to  known  and  computed17  (by  direct  numeri¬ 
cal  simulations)  properties  of  homogeneous  turbu¬ 
lence.  The  model  was  extended  by  Pope13  to  in¬ 
homogeneous  turbulence  and  validated  against 
several  self-similar  flows.  An  attempt  to  validate 
foe  model  for  a  piloted  jet  diffusion  flame  by  Nor¬ 
ris  and  Pope1®  was  not  very  successful  due  to  foe 
complexities  and  uncertainties  in  the  flow  studied. 
In  the  present  study,  foe  model  is  validated  against 
benchmark  quality  data  of  Takahashi  et  al20  for  a 
coaxial  swirling  jet  flow. 

A  unique  feature  of  the  laser  Doppler  ve- 
lodmeter  (LDV)  measurements  of  Takahashi  et  al2® 
is  that  conditional  velocity  statistics  are  measured 
by  seeding  individual  jets  one  at  a  time.  This  elim¬ 
inates  issues  of  velocity  bias  in  the  measurements. 
Conditional  turbulent  correlations  up  to  fourth  or¬ 
der  have  been  reported.  These  conditional  correla¬ 
tions  are  easily  computed  in  the  pdf  method,  and 
the  comparisons  against  data  serve  to  verify  that 
turbulent  transport  is  indeed  accurately  calculated 
in  the  pdf  method. 


A  brief  description  of  foe  method  pertinent  to 
foe  present  study  is  presented  here.  More  details 
can  be  obtained  from  previous  studies.2,12'16'1® 

The  equation  for  foe  evolution  of  foe  joint  pdf 
of  velocity,  dissipation,  and  other  scalars  is  mod¬ 
eled  and  solved  using  a  Lagrangian  viewpoint.  In 
the  pdf  transport  equation,  the  terms  to  be  modeled 
are  those  representing  the  effects  of  viscous 
dissipation,  the  fluctuating  pressure  gradient,  and 
molecular  diffusion.  All  other  processes  are  in 
closed  form  and  need  not  be  modeled.  The  solu¬ 
tion  is  performed  by  a  Monte  Carlo  algorithm.  In 
the  boundary-layer  algorithm  used  here,  the  joint 
pdf  at  each  step  in  foe  x-direction  is  represented  by 
a  large  number,  N,  of  notional  or  modeled  parti¬ 
cles.  At  foe  axial  position  x,  each  particle  has  the 
position  x*(x),  velocity  LJ*(x),  relaxation  rate  or  fre¬ 
quency  <o*(x),  and  a  scalar  value  c*(x).  The  dissipa¬ 
tion  model16,1®  is,  in  fact,  expressed  in  terms  of  to* 
(=  e  */  k )  rather  than  in  terms  of  the  particle  dissi¬ 
pation  e*  (foe  tilde  as  in  k,  denotes  ©-weighted 
means).  The  scalar  c*  is  a  passive  scalar  with  no 
source  terms  (i.e.,  c*  is  held  constant),  and  is  used 
to  tag  particles  by  their  jet  of  origin.  In  the  general 
step  from  x  to  x  +  Ax,  each  particle  evolves  over  a 
time  interval  At*  given  by: 

At*  =  Ax/U*i .  (1) 
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The  particle  velocity  and  relaxation  rate  evolve 
by  coupled  stochastic  equations  described  in  die 
next  section.  Means  and  other  correlations  are  ex¬ 
tracted  from  die  solution  by  forming  sums  of  par¬ 
ticle  properties  within  spatial  bins  and  fitting  cubic 
B-splines  to  these  sums. 

The  Stochastic  Velocity-Dissipation  Model 

The  details  on  die  construction  of  the  model 
and  die  choice  of  constants  are  discussed  exten¬ 
sively  by  Pope  and  Chen16  and  Pope.18  The  fol¬ 
lowing  notation  is  used  in  the  description  of  the 
models.  The  density  is  constant  in  the  present 
study  and  is  denoted  by  p.  The  Eulerian  instanta¬ 
neous  velocity  at  position  x  and  time  t,  12  (x,t)/  is 
written  in  terms  of  its  mean  <12  (x,t)>  and  fluctua¬ 
tion  u  (x,t)  as: 

12  =  <12>  +  u.  (2) 


entrained  into  a  turbulent  region  (<o»  >  0).  In  the 
present  study,  the  entire  flow  was  assumed  to  be 
turbulent  with  ou*  >  0  for  all  particles.  Neverthe¬ 
less,  h  «a>2  is  an  additional  source  term  for  the 
<fi»  equation  in  the  (non-Gaussian)  turbulent  re¬ 
gion  and  can  affect  the  results  if  Cgg  is  chosen  large 
enough. 

A  model  for  S©  that  is  consistent  with  the  k-e 
model  (and  the  model  for  <(o>  used  previ¬ 
ously1^!)  is:16 


S®j =  ■Qoi  Sjj  Sjj/  <(D>  +  C<n2/ 
where 


Sii 


1 

2 


9<Uj  >  |  d,<U)>" 
,  3xi  , 


(6) 


(7) 


According  to  the  model  for  inhomogeneous 
flows,  the  increment  to  to*  in  the  time  interval  dt  is 
given  by  (note:  angled  brackets  denote  means): 

d<o*  =  -<0*  <e»  (Sc8  +  C xCl)  dt  +  <©>2  h  dt 

+  c o *  (2C^  «d>  a2)1/2  dW,  (3) 

where 


Q  =  tn  tn  — 

V<fi>>,/  \<  00  >  <©>/ 


(4) 


and  dW  is  a  Weiner  process  with  properties  <dW> 
=  0  and  <dWdW>  =  dt.  The  constants  in  Equation 
3  have  been  assigned16  die  values  =  1.6  and  o2 
=  1.0.  According  to  the  model  (Equation  3)  the  Eu¬ 
lerian  mean  relaxation  rate  <co>  evolves  by 

d  <  €0  >  O 

— — —  =  -  <  to  >2  (Sj#  -  h).  (5) 

The  term  due  to  h  in  Equation  5  is  a  result  of 
die  extension  of  the  model  to  inhomogeneous  flows 
and  is  zero  in  homogeneous  turbulence.  The  term 
h  (described  in  detail  by  Pope18)  involves  a  con¬ 
stant  Cos  with  a  suggested  value  of  Cm3  =  1.0.  In 
die  present  study,  the  change  in  the  value  of  C©3 
over  a  range  of  0.0  to  5.0  produced  no  appreciable 
differences  in  the  results  obtained.  While  more  in¬ 
vestigation  is  needed  to  determine  the  magnitude 
of  Qo3,  the  term  due  to  h  played  an  insignificant 
role  in  the  present  calculations  for  the  values  of 
C<o3  used.  A  major  role  of  the  term  h  <e»2  is  to  in¬ 
crease  the  value  of  to*  for  a  nonturbulent  particle 
(of  =  0)  and  make  it  turbulent  (to*  >  0)  when  it  is 


It  should  be  noted  that  with  the  stochastic  ve¬ 
locity-dissipation  model,  there  is  no  need  to  restrict 
the  specification  of  S<o  in  terms  of  <co>  and  the 
mean  velocity  gradients.  Indeed,  as  more  insight  is 
gained  regarding  the  behavior  of  to*  (and  <co>),  Sw 
can  be  specified  to  include  other  correlations  ex¬ 
tracted  from  the  joint  pdf  of  to  and  velocity.  With 
the  specification  in  Equation  6,  the  resulting  equa¬ 
tion  for  <(0>  (Equation  5)  is: 

d<oo>  9<a>>  TT  3<to>  9<uf(o> 

- - - +<U;> - + - - - = 

dt  9t  3xj  9xj 

CuiSijSjj  -  C.2  <  to  >2  +h  <  to  >2.  (8) 

The  values  of  Qoi  and  Qo2  were  varied  over 
the  ranges  0.03  to  0.1  and  0.60  to  0.92,  respectively, 
in  the  present  study.  The  set  Cmi  =  0.04  and  C<n2  = 
0.70  was  found  to  result  in  the  best  agreement  with 
data  in  the  present  study.  These  values  are  still 
considered  to  be  tentative  until  the  model  is  vali¬ 
dated  against  several  more  flows.  The  values  C<oi 
=  0.04  and  C™?  =  0.09  were  used  by  Pope18  for  self¬ 
similar  flows.  In  that  study,  the  value  of  C®i  was 
chosen  to  yield  the  correct  value  of  the  von  Karman 
constant  (k  =  0.4)  in  the  constant-stress  wall  layer 
calculations,  and  the  value  of  C<o2  was  chosen  to  be 
consistent  with  the  k-e  model  and  not  varied. 

The  stochastic  model  for  the  velocity  increment 
over  a  time  interval  dt  is:18 

dU*  =  d^F>dt  +  Djdt  +  (CDk  to*)1/2  dWir  (9) 

1  p  9xi 
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where 


Di  =-^”+“'Co j<fi>>(^/k)ui  +  G^  uJ-jC0 
[  AJj1  (to  *  uj  -  <  mu j  >)  -  (k  /  k)A":  <  0)  >  uj],  (10) 

U:  =  Uj-<Ui>/  (ID 

and 

Aij  *  3<uj  up/<uk  uj(>.  (12) 

The  mean  pressure  is  <P>,  Di  is  the  drift  term, 

Ajj  is  the  normalized  Reynolds  stress  tensor  and 
dWi  is  an  isotropic  Weiner  process,  independent  of 
that  in  Equation  3.  In  general,  G*j  is  a  second-or¬ 
der  tensor  function  of  the  Reynolds  stresses  and 
mean  velocity  gradients.  The  choice  Gfj  =  0,  made 
here  and  in  References  18  and  20,  results  in  the 
evolution  equation*®  of  the  Reynolds  stresses 
(derived  from  Equation  9)  having  the  same  form  as 
the  Rotta's  model  in  conventional  Reynolds-stress 
closure.  For  the  Kolmogorov  constant  C0,  the 
value  Co  =  3.5  recommended  by  Pope  and  Chen*® 
is  used. 


Method  Based  on  Transport  Model  for  <co>  (<m>- 
Model) 


Before  the  development  of  the  stochastic  veloc¬ 
ity-dissipation  model,  some  previous  studies*2**® 
have  used  a  pdf  method  based  on  a  transport 
model  for  «o>,  henceforth  referred  to  as  the  <co> 
-model.  This  model  is  also  consistent  with  the  k-e 
model  and  can  be  derived  from  the  equations  for  k 
and  <e>  with  a  remodeling  of  the  turbulent 
convection  term.  The  equation  for  <©>  from  this 
model  is: 


3<to>  9<co>  _ _ 

-+  <  U|  >  — - 2Ch  (Ce1  - 1)  Sjj  Sjj  - 


dt 


dx, 


(Ce2-1)<<0>2  +  **  - 


°C 0 *4 


f  \ 

k  9<(o> 

<  (0  >  dx . 


,  (13) 


where  the  turbulent  convection  term  is  modeled  by 
gradient  diffusion: 


d  <  ut(0  > 

ax, 


cH  a  k  a<(o>> 

3x,  y<<n>  3x,  > 


(14) 


The  model  constants  are15  C^  =  0.09,  Cei  = 
156,  Cg2  =  1*88,  and  Oo  =  1.0.  The  evolution  equa¬ 
tions  for  «o>  horn  the  stochastic  model  (Equation 
8)  and  the  «o>-model  (Equation  13)  have  the  same 
form  except  that  h<o»2  is  an  additional  term 
(although  small  in  the  present  study)  in  Equation  8 
and  the  turbulent  convection  term  is  modeled  by 
Equation  14  in  Equation  13.  However,  the  constant 
2Cji(Qi  *  1)  =  0.1  is  significantly  different  from  the 
corresponding  constant  Qoi  =  0.04,  and  the  con¬ 
stant  (C{2  - 1)  =  0.88  is  different  from  Qb2  =  0.7. 

The  solutions  for  <co>  obtained  from  the  two  mod¬ 
els  are  significantly  different.  On  the  other  hand, 
the  stochastic  model  for  velocity  based  on  <co>  is 
also  significantly  different: 


9<P> 

- dt 

3x, 


*  3  • 

I+7C° 


\ 

y 


<  ©  >  u*  dt  +  (C„  k  <  (0  >)1/2  dW,. 


(15) 


The  differences  between  Equations  9  and  15  are 
that  the  random  term  (involving  dWj)  in  Equation 
9  contains  the  instantaneous  value  oo*  rather  than 
<(D>,  and  there  are  additional  drift  terms  in  Equa¬ 
tion  9  (see  Equation  10).  Further,  the  value  of  C0 
appropriate  for  this  model  (Equation  15),  deter¬ 
mined  by  Anand  and  Pope,22  is  C’0  =  2.1  (see  Pope 
and  Chen*®  for  further  discussion  on  these  differ¬ 
ences).  The  evolution  equation  for  <uiuj>  resulting 
from  Equation  15  also  corresponds  to  Rotta's  model 
in  the  conventional  Reynolds  stress  closure. 

In  the  present  study,  calculations  are  per¬ 
formed  using  both  the  stochastic  model  and  the 
<co>-model,  and  the  results  are  compared  against 
each  other  and  against  experimental  data. 

Algorithm  for  Swirling  Flows 

There  are  no  theoretical  limitations  for  comput¬ 
ing  swirling  flows  with  the  pdf  method.  Exten¬ 
sions  were  made  to  the  boundary-layer  algorithm 
to  be  able  to  compute  swirling  flows.  For  the  dis¬ 
cussions  to  follow,  the  notation  H  m  (Ui,  U2,  U3)  = 
(U,  V,  W)  s  (<U>  +  u,  <V>  +  v,  <W>  +  w)  is  used. 

In  addition,  due  to  the  geometry  of  the  flow  con¬ 
sidered,  polar  cylindrical  coordinates  (x,  r,  0)  are 
used.  The  main  extensions  to  the  algorithm  are  to 
calculate  and  include  the  axial  pressure  gradient 
and  to  include  all  terms  arising  due  to  <W>  *  0  in 
calculating  SjjSij  and  the  (radial  and)  axial  pressure 
gradients. 

The  mean  pressure  gradient  is  computed  by  in¬ 
tegrating  the  mean  radial  (or  lateral)  pressure  gra- 
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dient  across  the  radius  of  the  flow.  The  mean  ra¬ 
dial  pressure  gradient  is  given  by 


—  --  ---  =  —  (<  W  >2  +  <  w2  >  -  <  v2  >  - 
p  dr  r 

d<  v2  >. 

r-s->' 


(16) 


where  the  boundary  layer  assumptions  (neglect  of 
axial  gradients)  and  axisymmetry  (azimuthal  gra¬ 
dients  being  zero)  have  been  used  and  the  mean 
radial  velocity  and  its  gradients  are  neglected.  The 
mean  pressure,  as  a  function  of  r  at  a  given  x,  is 
given  by 


<  P(r ;  x)  >  =  <  P(Rq;x)  >  -  J 


r 


a<p> 

3r 


dr. 


(17) 


The  flow  considered  has  a  central  nonswirling 
jet  of  diameter  D  =  9.45  mm,  a  concentric  annular 
swirling  jet  of  diameter  Da  =  26.92  mm  surrounded 
by  a  low  velocity  coflow  in  a  test  section  of  150  mm 
diameter.  All  the  jets  are  of  air.  The  bulk-averaged 
axial  velocities  for  the  central  jet,  annulus,  and  co¬ 
flow  are  100, 20,  and  4  m/s,  respectively.  The 
Reynolds  number  for  the  central  jet  is  approxi¬ 
mately  21,400.  The  swirl  number  is  defined  as  the 
ratio  of  the  axial  flux  of  angular  momentum  to  that 
of  axial  momentum  divided  by  the  outer  radius: 

R. 

Jp<UxW>rdr 

S  =  -L- - .  (20) 

R0  |  p  <  U  >2  rdr 

R,  f 


where  Rq  is  the  radius  of  the  outer  boundary  of  the 
jet.  The  mean  axial  pressure  gradient  at  x  +  Ax,  as  a 
function  of  r,  is  then  deduced  by 


3<P(r;x  +  Ax)> 
3x 


<P(r;x  +  Ax)>-<P(r;x)> 
Ax 


.(18) 


Any  externally  imposed  axial  pressure  gradient 
can  be  incorporated  through  <P(R<>;  xj>  in  Equa¬ 
tion  17.  In  the  present  study,  <P(Rq;  x>  is  taken  to 
be  constant  (zero)  for  all  x.  The  axial  pressure 
gradient  is  important  in  calculating  the  evolution  of 
the  mean  axial  velocity  (especially  its  decay  on  the 
centerline)  when  appreciable  swirl  is  present  in  the 
flow. 


The  expression  for  Sij  Sjj  is: 
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Again,  the  terms  involving  <V>  in  Equation  19 
are  neglected  in  comparison  with  the  other  terms. 

Calculations 


Flow  Description 

Calculations  are  performed  for  a  constant- 
density  coaxial  swirling  jet  flow  for  which  bench¬ 
mark  quality  data  from  LDV  measurements  have 
been  reported  by  Takahaski  et  al.^° 


The  swirl  in  the  annular  jet  is  produced  by  a 
vaned  swirler  with  a  vane  helix  angle  of  30  deg,  lo¬ 
cated  96  mm  upstream  of  the  nozzle  exit.  The 
swirl  number  for  the  annular  jet,  Sa,  based  on  Rj  = 
R  =  D/2  and  Rq  *  Ra  =  Da/2,  calculated  from  mea¬ 
sured  data  15  mm  downstream  of  the  nozzles,  is 
0.41.  However,  the  overall  swirl  number,  S0v/ 
based  on  the  central  and  annular  jets  (with  Rj  =  0 
and  Rq  =  Ra)  is  0.09  due  to  the  large  bulk  axial  ve¬ 
locity  of  the  central  jet.  Since  the  overall  swirl 
number  is  low,  the  inclusion  of  the  axial  pressure 
gradient  did  not  have  a  significant  effect  on  the  re¬ 
sults.  However,  with  another  test  problem^  with 
an  overall  swirl  number  of  approximately  0.4,  the 
inclusion  of  the  axial  pressure  gradient  was  neces¬ 
sary  to  obtain  the  correct  (measured)  spreading  of 
the  mean  axial  velocity  profile. 

Initial  Conditions 

Initial  conditions  for  the  computations  are  pre¬ 
scribed  from  experimental  data.  The  first  mea¬ 
surement  station  x  *  15  mm  (x/D  =  0.16)  is  taken 
as  the  initial  plane  for  the  calculations.  The  initial 
velocity  pdf  is  prescribed  to  be  joint  normal  with 
the  mean  and  covariances  taken  from  linearly  in¬ 
terpolated  experimental  data.  The  initial  pdf  of  re¬ 
laxation  rate  for  the  stochastic  model  is  taken  to  be 
log-normal  [ln(o>*/«*»)  has  a  normal  distribution 

with  mean  and  variance  —Jo2  and  a2,  respec¬ 
tively],  in  accordance  with  the  construction  of  the 
model.  The  specification  of  the  initial  <o»  profile 
for  the  stochastic  model  is  discussed  later  in  this 
section. 

The  initial  profile  of  <a»  for  the  <a»-model  is 
derived  from  experimental  data  using  the  expres¬ 
sion  based  on  the  k-e  model: 
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where  equivalent  velocity  gradient  and  shear  stress 
are  given  by: 

fa<u>l 


ra<u>>|2  p<w> 

i.  3r  j  l  dr 


<W> 


|l/2 


,  (22) 


tions)  was  poorer.  In  particular  the  turbulent  ki¬ 
netic  energy  was  lower  than  in  the  experiments  in 
this  region.  This  would  have  an  effect  on  the  re¬ 
sults  downstream,  especially  the  spreading  rate. 
The  results  suggested  that  the  appropriate  value  of 
inlet  «a>  for  the  stochastic  model  is  about  70%  of 
the  value  used  for  the  «o>  model,  keeping  in  mind 
that  the  stochastic  model  has  fluctuations  in  co  and 
noting  the  differences  in  the  velocity  models  dis¬ 
cussed  previously.  Hence  the  inlet  <e»-f6r  the 
stochastic  model  is  specified  to  be  0.7  times  the 
value  obtained  from  Equation  24.  Again,  the  alter¬ 
native  specification  of  0.7  times  the  value  obtained 
from  Equation  24  produced  nearly  the  same  re¬ 
sults. 


and 


<uv>eq  =  (<uv>2  +  <uw>^)l/2. 


(23) 


The  initial  «d>  profile  was  also  specified  from 
Equation  21  with  the  assumption  of  local  equilib¬ 
rium  of  turbulence  energy  in  the  shear  layer  that 
yields  <uv>eq  =  k  so  that: 


(24) 


Although  Equations  21  and  24  result  in  slightly 
differing  profiles  of  initial  «o>,  both  specifications 
produced  nearly  the  same  results  at  the  down¬ 
stream  stations.  For  the  results  with  the  «o>- 
model  reported  here,  the  specification  in  Equation 
24  was  used. 


The  specification  of  <to>  for  the  stochastic 
model  needs  some  clarification.  The  production  of 
<to>  in  the  stochastic  model,  controlled  by  C^i  = 
(0.04)  (Equation  6),  is  only  40%  of  the  production 
with  the  <co>  model  [Equation  13, 2C^  (Qi-1)  = 
0.1].  The  dissipation  of  <w>  is  also  lower  in  the 
stochastic  model,  but  is  about  80%  of  that  with  the 
«D>— model  [C(o2  =  0.70  versus  (Ce2  -1)  =  0.88]. 
Hence,  with  the  two  models  producing  nearly  the 
same  velocity  statistics,  the  stochastic  model  will 
yield  lower  values  of  <o»,  especially  in  the  near¬ 
nozzle  region  where  production  is  dominant  due  to 
large  velocity  gradients.  This  fact  is  borne  out  by 
the  results  presented  in  the  next  section.  With  the 
specification  of  the  same  inlet  <co>  for  both  the 
models,  the  stochastic  model  yielded  lower  <to> 

(by  about  30%)  than  the  «b>-  model  at  the  down¬ 
stream  locations  (x/D  =  0.5  to  10  approximately). 
The  agreement  with  data  of  the  turbulent  correla¬ 
tions  dose  to  the  nozzle  (x/D  =  1.06  and  2.65  sta- 


To  utilize  the  opportunity  provided  by  the 
stochastic  model  to  specify  additional  initial  condi¬ 
tions  based  on  joint  velodty-relaxation  rate  statis¬ 
tics,  initial  conditions  on  <vco>  were  specified  such 
that  the  correlation  coeffident  for  <v©>  is  the  nega¬ 
tive  (from  wall  boundary-layer  considerations)  of 
that  of  <uv>.  Although,  this  specification  did  not 
have  a  major  effect  on  the  results,  it  underscores 
the  fact  that  more  information  is  contained  in  and 
can  be  supplied  to  the  stochastic  model. 

Conditional  Statistics 

As  mentioned  previously,  LDV  measurements 
for  the  flow  were  made  by  seeding  individual  jets 
one  at  a  time.  The  complete  flow  field  would  be 
mapped  by  seeding  a  particular  jet  and  the  proce¬ 
dure  would  be  repeated  for  the  remaining  jets.  The 
experiments  were  shown  to  be  repeatable.  This 
technique  measures  the  conditional  pdfs  (i.e.,  con¬ 
ditional  upon  the  jet  of  origin)  and  avoids  errors 
due  to  velocity  bias  that  will  be  significant  due  to 
the  large  differences  in  the  bulk  velocity  among  the 
jets.  However,  the  unconditional  pdf  cannot  be 
deduced  from  these  measurements  without  the 
additional  measurement  of  the  fraction  of  time  the 
fluid  from  each  jet  spends  within  the  sample  vol¬ 
ume.  The  advantage  with  the  pdf  method  is  that 
the  conditional  pdfs  are  easily  calculated  in  addi¬ 
tion  to  the  unconditional  pdf.  Conventional  turbu¬ 
lence  models  can  only  calculate  (a  limited  number 
of)  unconditional  correlations  and  cannot  calculate 
conditional  quantities  without  requiring  additional 
modeling. 

The  comparison  between  conditional  mea¬ 
surements  and  calculations  sheds  valuable  light  on 
the  detailed  performance  on  the  model  and  on  the 
dynamics  of  individual  jets.  In  the  present  study, 
good  agreement  between  calculated  and  measured 
conditional  quantities  also  implies  that  the  calcu- 
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la  ted  unconditional  quantities  are  correct  since  no 
additional  modeling  is  required  to  derive  the  con¬ 
ditional  quantities. 

In  the  calculations,  the  particles  are  tagged  at 
the  initial  plane  by  assigning  the  value  of  c*  to  be  1, 
2,  or  3  depending  on  their  jet  of  origin.  The  value 
of  c*  does  not  change  throughout  the  calculation. 
Conditional  quantities  are  computed  for  output  by 
forming  sums  over  particles  with  the  same  value  of 
c*  within  a  spatial  cell  and  fitting  B-splines  to  these 
sums. 

Computational  Parameters 

The  total  number  of  particles  used  in  the  simu¬ 
lations  is  150,000.  A  large  number  is  used  to  re¬ 
duce  statistical  error  especially  in  the  calculation  of 
the  third  and  fourth  order  moments  (mean  and 
second  order  quantities  can  be  calculated  to  within 
a  few  percent  accuracy  with  as  low  as  25,000  parti¬ 
cles).  The  number  of  cells  used  for  forming  sums  is 
44  and  the  number  of  basis  functions  used  for  13- 
splines  is  20. 


x/P  =  1 .06  5.29 


<U>/<U>oc 


The  stochastic  model  calculations  took  approx¬ 
imately  38  minutes  of  CPU  time  to  march  from  x/D 
=  1.06  to  x/D  =  26.5  (in  about  540  steps)  on  a  CRAY 
YMP.  The  <co>-model  took  about  28  minutes  and 
470  steps.  No  significant  effort  has  been  made  to 
optimize  the  computer  code  for  speed.  The  storage 
required  was  about  7  megawords. 


Figure  1.  Radial  profiles  of  conditional  mean  axial 
velocity  from  pdf  (stochastic  velocity-dissipation 
model)  calculations  compared  against  data  at  vari¬ 
ous  axial  stations.  Lines  denote  calculations  and 
symbols  denote  experimental  data:  solid  line,  CD- 
central  jet  fluid;  short  dashed  line,  "-annular  jet 
fluid;  long  dashed  line,  V-coflow  fluid. 


Results  and  Discussion 

The  measured  mean  axial  velocity  on  the  cen¬ 
terline  of  the  jet  at  the  x  =  1.5  mm  station,  <U>oc 
(=129.7  m/s),  is  used  for  normalizing  all  the  veloc¬ 
ity  statistics  presented.  The  radial  distance  is  nor¬ 
malized  by  the  radius  of  the  central  jet  (R)  and  the 
axial  distance  is  normalized  by  the  central  jet  di¬ 
ameter  (D). 

Measurements20  were  taken  across  the  diame¬ 
ter  of  the  flow.  However,  in  the  figures  presented 
here,  the  data  on  either  side  of  the  centerline  are 
shown  on  the  same  side.  This  gives  an  indication 
of  the  symmetry  of  the  flow  and  the  amount  of 
scatter  in  the  data.  Measurements  have  been  re¬ 
ported  for  x/D  =  0.16  (initial  plane),  1.06, 2.65, 5.29, 
7.94, 15.9,  and  265.  Figure  1  shows  the  conditional 
mean  axial  velocities  from  the  stochastic  model  at 
four  different  axial  stations.  The  calculations  are  in 
very  good  agreement  with  data.  The  conditional 
velocities  show  the  momentum  transfer  between 
the  jets  as  they  mix.  The  annular  jet  axial  velocity 
is  accelerated  as  the  annular  fluid  moves  towards 


the  centerline  and  the  central  jet  decelerates.  Simi¬ 
larly,  the  coflow  fluid  accelerates  as  it  moves  in. 
The  conditional  velocities  differ  significantly  from 
each  other  in  the  developing  region,  and  these  dif¬ 
ferences  become  insignificant  far  downstream  (x/D 
=  26.5). 


The  corresponding  mean  axial  velocity  profiles 
from  the  <(o>-model  are  shown  in  Figure  2.  The 


tional  velocities  show  the  same  relative  trends  and 
the  agreement  with  data  over  most  of  the  radius  is 
quite  good.  The  spreading  rate  (and  the  centerline 
velocity  decay)  can  be  adjusted  by  changing  the 
values  of  Qi  and  Q2,.  However,  tests  indicated 
that  changing  the  constants  to  get  a  better  agree¬ 
ment  on  the  spreading  rate  would  also  lower  the 
turbulent  kinetic  energy  (i.e.,  <u2>,  <v2>,  and 
<w2>)  and  results  in  a  poorer  agreement  with  data 
for  those  and  other  turbulent  correlations. 


profiles  from  the  <co>-model  show  greater  decay 
the  centerline  velocity  than  those  shown  by  the 
stochastic  model  and  data.  However,  the  condi- 
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Figure  2.  Radial  profiles  of  conditional  mean  axial 
velocity  from  pdf  (<G»-model)  calculations  com¬ 
pared  against  data  at  various  axial  stations  (see 
Figure  1  for  legend). 

The  profiles  of  <(o>  obtained  from  the  stochas¬ 
tic  and  «o>-model  are  shown  in  Figure  3.  As  ex¬ 
plained  earlier,  the  <o»  profiles,  especially  the 
peak  values,  from  the  stochastic  model  are  consid¬ 
erably  lower  than  those  from  the  <a»-model  in  the 
developing  region  (x/D  <  7.94).  In  this  region  the 
production  of  «b>  plays  a  major  role  due  to  the 
large  mean  velocity  gradients.  The  production  is 
lower  in  the  stochastic  model  due  to  the  choice  of 
the  constants.  Further  downstream,  production  is 
less  and  the  dissipation  of  <o»  plays  a  significant 
role.  Since  the  dissipation  is  lower  in  the  stochastic 
model  than  in  the  <e»-model  it  can  be  seen  that 
the  <to>  values  are  higher  for  the  stochastic  model 
at  x/D  =  26.5.  In  fact,  the  two  profiles  crossover  at 
about  x/D  =  16. 

The  profiles  of  the  mean  swirl  velocity  are 
shown  in  Figure  4.  The  decay  of  the  swirl  velocity 
is  well  predicted  at  all  stations  including  the  down¬ 
stream  stations  not  shown  here.  The  increase  in 
swirl  velocity  of  the  central  jet  and  coflow  fluids  as 
they  move  towards  the  annular  swirling  jet  is 
clearly  depicted  by  both  the  measured  and  calcu¬ 
lated  profiles.  The  profiles  of  <W>  calculated  by 
the  «a>-model  are  also  in  good  agreement  with 
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Figure  3.  Profiles  of  mean  relaxation  rate  <co>  from 
stochastic  model  calculations  (solid  lines)  com- 
;  pared  with  those  from  <co>-model  calculations 
(dashed  lines)  at  various  axial  stations. 

data.  The  negative  values  of  <W>  shown  by  data 
at  x/D  =  5.29  and  7.94  are  nonphysical  and  indica¬ 
tive  of  some  experimental  error. 

Since  most  significant  differences  among  the 
conditional  quantities  exist  in  the  near  nozzle  re¬ 
gion,  x/D  <  7.94,  only  selected  stations  in  this  re¬ 
gion  are  shown,  in  the  interest  of  clarity,  in  the  fig¬ 
ures  to  follow.  The  results  are  in  good  agreement 
with  data  at  all  stations.  The  profiles  of  turbulent 
kinetic  energy  from  the  stochastic  model  are  shown 
in  Figure  5  for  x/D  =  5.29  and  7.94. 

Figure  6  shows  the  unconditional  and  condi¬ 
tional  mean  axial  velocity  and  turbulent  kinetic  en¬ 
ergy  at  x/D  =  5.29.  The  unconditional  mean  veloc¬ 
ity  lies  within  the  conditional  velocities  at  any 
given  radius  since: 
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Figure  4.  Profiles  of  conditional  mean  swirl  veloc¬ 
ity  <W>  from  stochastic  model  calculations  com¬ 
pared  against  data  (see  Figure  1  for  legend). 
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x/D  =  5.29  7.94 


IOOk/<U>2oc 

Hgure5.  Profiles  of  conditional  turbulent  kinetic 
energy  from  stochastic  model  calculations  com¬ 
pared  against  data  (see  Figure  1  for  legend). 
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Figure  6.  Profiles  of  conditional  and  unconditional 
mean  axial  velocity  and  turbulent  kinetic  energy  at 
x/D  =  5.29  from  stochastic  model  calculations. 
Solid  line — -jet  fluid;  short  dashed  line — annular  jet 
fluid;  long  dashed  line — coflow;  extra  long  dashed 
line— unconditional. 

<U>  =  Kj  <U>j  +  Ya  <U>a  +  Yc  <U>C  (25) 

where  file  subscripts  j,  a,  and  c  denote  the  central 
jet,  annular  jet,  ami  coflow,  respectively,  and  y  is 
file  fraction  of  file  time  that  the  fluid  from  the 
particular  jet  is  found  at  the  measurement  location 
(note  that  7s  are  less  than  unity  and  add  up  to 
unity).  However,  the  unconditional  turbulent 
correlations,  in  general,  need  not  lie  within  the 
conditional  quantities.  Expressions  similar  to 
Equation  25  can  be  derived  for  each  of  the 
correlations  that  will  have,  in  addition  to  a  linear 
combination  of  the  conditional  correlations,  terms 
involving  differences  between  the  conditional 
mean  velocities  and  other  conditional  quantities.  It 
is  easier  to  visualize  that  while  the  conditional 
fluctuations  are  calculated  from  the  corresponding 
conditional  mean,  the  unconditional  fluctuations 
are  with  respect  to  the  unconditional  mean.  For 


example,  in  a  region  where  the  conditional 
velocities  are  significantly  different,  the  root  mean 
square  (rms)  fluctuation  about  the  unconditional 
mean  can  be  significantly  larger  than  the  rms 
fluctuations  about  the  respective  conditional 
means.  This  fact  is  illustrated  by  the  kinetic  energy 
profiles  in  Figure  6. 

The  kinetic  energy  profiles  from  the  <co>- 
model  (not  shown  here)  are  also  in  good  agreement 
with  data,  but  some  differences  in  the  details  exist 
between  the  results  from  the  two  models  as  illus¬ 
trated  by  the  next  two  figures.  Figures  7  and  8 
show  the  conditional  axial  velocity  variance  from 
the  stochastic  and  <o»-models  respectively.  Both 
models  produce  results  that  have  the  same  level  of 
agreement  with  data.  However,  the  shift  in  the  lo¬ 
cation  of  the  peaks  for  the  central  and  annular  jets 
is  better  captured  by  the  stochastic  model  than  the 
<o»-model.  This  may  be  due  to  the  fact  that  the 
turbulent  transport  of  to  is  modeled  (by  gradient 
diffusion)  in  the  <(0>-model  while  it  need  not  be 
modeled  with  the  stochastic  model. 


x/D  =  2.65  5.29 
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Figure  7.  Profiles  of  conditional  axial  velocity  vari¬ 
ance  from  stochastic  model  calculations  compared 
against  data  (see  Figure  1  for  legend). 
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Figure  8.  Profiles  of  conditional  axial  velocity  vari¬ 
ance  from  <oo>-model  calculations  compared 
against  data  (see  Figure  1  for  legend). 
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This  discrepancy  in  the  <u2>  profile  is  further 
amplified  in  the  <u2v>  profiles  shown  in  Figures  9 
and  10.  The  quantity  <u2v>  can  be  interpreted  as 
the  radial  turbulent  flux  of  <u2>.  The  quantity 
<u*v>  will  have  opposite  signs  on  either  side  of  the 
location  (approximately)  of  the  peak  of  <u2>  since 
turbulence  will  disperse  the  peak  in  opposite  direc¬ 
tions  from  the  peak.  Figure  9  shows  that  the 
stochastic  model  shows  good  agreement  with  data 
while  the  <o»-model  calculations  are  not  in  as 
good  an  agreement  (Figure  10). 


In  Reynolds-stress  closures,  all  third-order  cor¬ 
relations  (such  as  <u2v> )  are  modeled,  usually  by 
gradient  diffusion.  In  the  present  study,  the  flow 
under  consideration  was  calculated  using  the  RSM 
of  Daly  and  Harlow2*  that  is  widely  used.  The 
details  of  the  code  and  models  used  are  presented 
in  Reference  25.  The  model  for  <u2v>  is: 


One  can  draw  an  analogy  between  <u2v>  (the 
transport  of  the  velocity  variance  <u2>  in  the  v-di- 
rechon)  and  <^v>  (the  transport  of  «|>2>,  variance 
of  a  scalar  ®  in  the  v-direction).  It  is  important  to 
calculate  scalar  transport  (both  mean  and  fluctua¬ 
tion  of  the  scalar)  accurately,  especially  in  reacting 
flows,  since  it  significantly  affects  local  heat  release 
rates.  It  is  even  more  critical  when  there  is  mixing 
between  different  streams  (such  as  fuel  and  oxi¬ 
dizer  in  diffusion  flames). 


Figure  9.  Profiles  of  (conditional)  triple  correlation 
<u2v>  from  stochastic  model  calculations  com¬ 
pared  against  data  (see  Figure  1  for  legend). 
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Figure  10.  Profiles  of  (conditional)  triple  correla¬ 
tion  <u2v>  from  <co>-model  calculations  compared 
against  data  (see  Figure  1  for  legend). 


where  Cg  =  0.22.  The  profile  of  (unconditional) 
<u2v>  calculated  from  Equation  26  with  the  right- 
hand  side  quantities  taken  from  the  Reynolds- 
stress  calculations  (short  dashed  line)  is  compared 
against  the  unconditional  <u2v>  profile  from  the 
stochastic  model  calculations  (solid  line)  in  Figure 
11.  The  comparison  shows  that  the  RSM  calcula¬ 
tions  severely  underpredict  the  transport  of  <u2> 
especially  at  the  stations  close  to  the  nozzle.  The 
magnitude  of  <u2v>  is  lower  because  the  RSM  cal¬ 
culations  also  grossly  underpredict  <u2>  in  this  re¬ 
gion.  In  addition  the  <u2v>  profiles  from  the  RSM 
calculations  change  sign  at  a  different  radial  loca¬ 
tion  than  the  pdf  calculations  since  the  peaks  of 
<u2>  are  also  shifted  in  the  RSM  calculations. 

These  differences  get  larger  further  downstream. 

A  better  test  of  the  gradient  diffusion  hypothesis 
would  be  to  use  the  pdf  results  for  the  right-hand 
side  in  Equation  26.  These  profiles  (long  dashed 
lines  in  Figure  11)  also  show  differences  both  in 
magnitude  and  the  location  where  the  profiles 
change  sign.  Though  these  discrepancies  are  not 
dramatic  for  the  flow  considered  here,  they  under¬ 
score  the  fact  that  the  gradient  diffusion  hypothesis 


x/D -1.06  2.65  5.29 


Figure  11.  Profiles  of  (unconditional)  triple  corre¬ 
lation  <u2v>  from  stochastic  model  calculations 
(solid  line);  from  Reynolds-stress  model  (Equation 
26)  using  results  from  Reynolds-stress  closure  cal¬ 
culations  (short  dashed  line);  and  from  model 
(Equation  26)  using  results  from  stochastic  model 
calculations  (long  dashed  line). 
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could  lead  to  different  results  and  in  the  raso  of  re¬ 
acting  flows  where  counter-gradient  diffusion  has 
been  demonstrated2^2®  the  results  would  be  erro¬ 
neous. 

It  should  be  noted  that  even  with  the  sralar  pdf 
method2*26  in  which  the  joint  pdf  of  scalars  only  is 
considered,  the  turbulent  flux  of  the  scalar  pdf  is 
modeled  by  gradient-diffusion.  While  the  scalar- 
pdf  method  is  a  useful  tool  for  the  development  of 
reduced  reaction  mechanisms,  since  reaction  still 
appears  in  closed  form  in  the  method,  the  full  po¬ 
tential  of  the  pdf  method  can  only  be  realized  by 
tite  velocity-scalar  joint  pdf  method. 

Figure  12  shows  the  conditional  <u®>  profiles 
at  x/D  =  2.65  and  529.  The  dramatic  differences 
between  the  conditional  quantities  is  well  predicted 
by  the  calculations.  The  results  from  the  «o>- 
model  show  similar  discrepancies  as  seen  for 
<u2v>.  Figure  13  shows  the  (conditional)  fourth 
order  correlation  <w4>  from  the  stochastic  model 
at  x/D  =  529  and  7.94.  Considering  that  fourth  or- 


X/D  =  2.65  5.29 


-1.20  -.80  -.40  0.0  .40  .80  -1 .20 -.80  -.40  0.0  .40  .80 


1000<u3>/<U>2oc 

Figure  12.  Profiles  of  conditional  <u®>  from 
stochastic  model  compared  against  data  (see 
Figure  1  for  legend). 

x/D  =  5.29  7.94 


3. 


1 0000<w4>/<U>4oc 


Figure  13.  Profiles  of  (conditional)  fourth-order 
correlation  <w4>  from  stochastic  model  calcula¬ 
tions  compared  against  data  (see  Figure  1  for 
legend). 


der  conditional  statistics  are  being  measured  and 
computed,  the  results  are  in  very  good  agreement 
with  data.  The  corresponding  profiles  for  the  <oo>- 
model  are  shown  in  Figure  14.  These  profiles  are 
also  in  good  agreement  with  data.  However  the 
radial  separation  between  the  peak  locations  for  the 
central  jet  and  annulus  at  x/D  =  5.29  are  better 
predicted  by  the  stochastic  model. 

Conclusions 

Swirling  flow  calculations  have  been  made 
with  the  pdf  method.  The  newly  developed 
stochastic  velocity-dissipation  model  has  been 
validated  against  benchmark  quality  data  fro  a 
coaxial  swirling  jet  flow.  Detailed  comparisons  be¬ 
tween  measured  and  calculated  conditional  means 
and  correlations  up  to  fourth  order  have  been 
made.  The  calculations  are  in  very  good  agreement 
with  data.  The  present  study  suggests  the  values 
Qol  =  0.04  and  C<q2  =  0.7  for  the  two  constants  that 
significantly  affect  the  evolution  of  the  mean  relax¬ 
ation  rate  <co>. 

Swirling  flows  were  calculated  for  the  first  time 
using  the  pdf  method.  While  extensions  to  the  al¬ 
gorithm  were  made  to  include  the  axial  pressure 
gradient  and  additional  terms  due  to  nonzero 
mean  swirl  velocity,  no  specific  modifications  were 
made  to  the  <o*  -  equation  or  elsewhere  on  account 
of  swirl. 

Calculations  were  also  made  with  the  pdf 
method  using  a  transport  model  for  <co>.  These  re¬ 
sults  are  also  in  reasonably  good  agreement  with 
data,  although  the  stochastic  model  better  predicts 
some  of  the  details  of  turbulent  transport. 

The  calculated  results  from  the  stochastic 
model  were  compared  against  gradient-diffusion 


Figure  14.  Profiles  of  (conditional)  fourth  order 
correlation  <w^>  from  «o>-model  compared 
against  data  (see  Figure  1  for  legend). 
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models  used  by  conventional  Reynolds-stress  clo¬ 
sures,  and  the  latter  were  found  to  be  deficient.  By 
extension,  the  scalar  pdf  method  is  also  subject  to 
the  same  limitations  since  gradient  diffusion  mod¬ 
eling  is  used  for  turbulent  transport  of  the  scalar 
pdf. 

The  pdf  model  used,  with  G(j  =  0  is  not  very 
sophisticated  in  comparison  with  available 
Reynolds-stress  closures.  It  corresponds  to  Rotta's 
model  with  no  rapid  pressure  terms.  It  is  remark¬ 
able,  therefore,  that  this  level  of  agreement  with 
data  is  achieved  even  with  this  simple  choice. 
Clearly,  a  more  general  model  for  G  S  should  be 
investigated  in  the  future. 

The  stochastic  velocity-dissipation  model  pro¬ 
vides  complete  closure  for  the  joint  velocity-scalar 
pdf  method  and  eliminates  the  need  to  supply  a 
time  scale  for  the  solution  of  the  pdf  transport 
equation.  The  stochastic  model  also  contains  more 
information  than  that  based  on  <o»,  and  allows  the 
flexibility  to  build  in  more  information  than  cur¬ 
rently  used  in  further  developing  the  model. 
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